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Abstract. We study neutralino sectors in extensions of the MSSM that dynamically generate the /l-term. The extra neutralino 
states are superpartners of the Higgs singlets and/or additional gauge bosons. The extended models may have distinct lightest 
neutralino properties which can have important influences on their phenomenology. We consider constraints on the lightest 
neutralino from LEP, Tevatron, and (g — 2)^ measurements and the relic density of the dark matter. The lightest neutralino 
can be extremely light and/or dominated by its singlino component, which does not couple directly to SM particles except 
Higgs doublets. 
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INTRODUCTION 

Extensive studies of the Minimal version of the Super- 
symmetric SM (MSSM) show that its lightest neutralino 
has the right ranges of mass and interaction strength to 
be a good cold dark matter (CDM) candidate [2]. Al- 
though the MSSM may be the optimal low energy Su- 
persymmetric model with minimal extension of the fields 
and symmetry, the model may not fully describe the TeV 
scale physics. It is important to see if other versions of 
the Supersymmetric SM can also give acceptable CDM. 

We consider various extended MSSM models that re- 
solve the jj. -problem [3] and compare their lightest neu- 
tralino properties to that of the MSSM. In these beyond- 
MSSM models, at least one Higgs singlet is commonly 
present and generates an effective jti parameter when the 
associated symmetry is broken at the EW or TeV scale. 
Because of the superpartners of the Higgs singlets or the 
extra gauge boson, the neutralino sector in these models 
may be significantly different from that of the MSSM. 
We investigate the mass and the coupling of the lightest 
neutralino of each model allowed by the model parame- 
ters and the experimental data. 

MODELS 

The extended MSSM models that we consider are 
the Next-to-Minimal Supersymmetric SM (NMSSM) 
[4], the Minimal Non-minimal Supersymmetric SM 
(MNSSM) a.k.a. the nearly Minimal Supersymmetric 
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SM (nMSSM) [5], the £7(1 /-extended Minimal Super- 
symmetric SM (UMSSM) [6], and the /-extended 
Supersymmetric SM with a secluded t/(l)'-breaking 
sector (S-model) [7]. All of these extended models pre- 
vent the /I -term (/1H1//2) and allow an effective jti-term 
(SH1H2) through a vacuum expectation value (VEV) (S) 
of a Higgs singlet associated with a new symmetry. 

The NMSSM assumes a discrete symmetry Z3 to 
avoid the /x-term, but allows an S 3 term in the superpo- 
tential. The NMSSM is one of the simplest extensions of 
the MSSM, but the Z3 symmetry predicts domain walls 
which are not observed [8]. The nMSSM was devised to 
avoid the domain wall problem while maintaining a dis- 
crete symmetry. The term aS in the nMSSM is a loop- 
generated tadpole term that breaks the discrete symme- 
try. The UMSSM uses an Abelian gauge symmetry in- 
stead of a discrete symmetry and thus is free from the do- 
main wall problem. The new gauge symmetry introduces 
a new gauge boson Z'. The S-model was introduced to 
resolve tension between the EW scale /i e ff and the heavy 
Z' (up to multi-TeV scale). It is basically the extension 
of the UMSSM with 3 additional Higgs singlets to pro- 
vide additional contributions to the Z' mass while keep- 
ing /i e ff = h s (S) at the EW scale. 

The models are listed in Table 1 with adopted sym- 
metries, superpotentials of Higgses, and the Higgses and 
neutralinos. The other parts of the superpotentials are the 
Yukawa terms of the MSSM and possible extra terms 
related to the exotic chiral fields in the 1/(1)' mod- 
els needed to cancel anomalies. In the /-extended 
model, the addition of one Higgs singlet does not give an 
additional CP odd Higgs since a goldstone boson is ab- 
sorbed to be the longitudinal mode of the massive f/(l/ 
gauge boson, Z'. 

For easy comparisons, we use the common notation of 



TABLE 1. Higgses and Neutralinos of the MSSM and its extensions 
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h s for the coefficient of the SH\H 2 term in each model. 
In every model the dynamically generated effective fi 
parameter is given by 



H eS = h s (S) 



(1) 



and therefore the VEV of the Higgs singlet or the sym- 
metry breaking scale needs to be at the EW/TeV scale. 



NEUTRALINO MASS MATRICES 

With the superpartners of the Higgs singlet (S) and the 
U(l)' gauge boson (Z'), the UMSSM has 6 neutralinos 
(MSSM components + S + Z') and its mass matrix (M^o) 
in the basis of {§, W 3 , H®, H$, S, Z'} is given by 
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and (//?) = 



gz'Q'^ (S) = ^ a„ u vu i = -j. 
v ~ 246 GeV. The gauge couplings are gi 



with 



v l + v 2 = 

e/cosOw and g 2 = e/sin0iy. For the U(l)' gauge cou- 
pling constant g z i and the U(l)' charge Q' , we use the 
Grand Unification Theory (GUT) motivated gauge cou- 
pling and the 77 -model charge assignments in our numer- 
ical analysis [1]. The U(l)' charges should satisfy 

Q' Hi +Q'h 2 ^° Q'hi +Q'h 2 + Q's = Q (2) 

in order to replace the /i-term with the effective jti-term 
dynamically generated by the Higgs singlet S. 

The first 5x5 submatrix corresponds to the nMSSM 
limit (or NMSSM limit in the special case K = 0) that can 
be obtained by taking M v > ^(EW). With a (5,5) entry 
of \flKs from S 3 term in the superpotential, this would 
be the NMSSM limit. 

The first 4x4 submatrix corresponds to the MSSM 
limit, which can be obtained by taking 5 ^> <^(EW) on 
top of the nMSSM limit. 

The S-model has 9 neutralinos (UMSSM compo- 
nents + §1, S 2 , S3), and its mass matrix has 3 more 



columns/rows added to UMSSM neutralino mass matrix. 
The UMSSM limit can be realized by taking 51,2,3 3* 
0(EW) with X s comparable to gauge couplings. How- 
ever, the most realistic case is for small X s and large 51,2,3 
[7], in which case four of the neutralinos, consisting al- 
most entirely of Z', S\, S 2 , and §3, essentially decouple 
from the other states. Since the full 9x9 matrix has a 
number of free parameters, we consider only this decou- 
pling limit when we discuss the light neutralinos, where 
there are 5 neutralinos with masses and compositions the 
same as the nMSSM 2 . 



DIRECT CONSTRAINTS 

The diagonalization of the neutralino mass matrix is 
accomplished via a unitary matrix N as 



N T M x0 N = Diag(M z o,M x o, • • •)• 



(3) 



The singlino (5) composition of the lightest neutralino 
{%]) 1S l-^isl 2 - Here we evaluate the bounds on the light- 
est neutralino mass (M o) and the singlino component 

(|Ms| 2 ) in the MSSM, the NMSSM, the nMSSM (also 
the decoupling limit of the S-model), and the UMSSM 
with gaugino mass unification. 

We require that the tree-level masses satisfy the direct 
LEP limits of M v ± > 104 GeV, and r„ n < 2.3 MeV 
(95% C.L.). Bounds from naturalness and perturbativ- 
ity constraints [7, 10, 11] are also imposed on the cou- 
plings of 0.1 < h s < 0.75 and sjh 1 + K 2 < 0.75 (for the 
NMSSM) 3 . The LEP2 Higgs mass bound of m h > 114 
GeV does not apply directly to the extended MSSM mod- 



2 The four decoupled neutralinos typically consist of one heavy pair 
involving the Z' and one linear combination of Si , 52 , S3 , as well as two 
more states associated with the orthogonal combinations of Si^.3 [7, 9]. 
The latter can be light, or even be the lightest neutralino in limiting 
cases. We do not consider that possibility here. 

3 An exact h s bound (and its source) may be a little different depending 
on models, but we assume a common bound for easy comparison. A 
lower bound on the NMSSM | k\ is fuzzy except that K ^ should 
be satisfied to avoid an unacceptable Peccei-Quinn symmetry. We set 
I If I > 0.1 as a lower bound; results for a smaller \k\ can be described 
by an interpolation of the nMSSM result (K —> limit). 



(a) Singlino compositions after direct constraints (b) Allowed mass ranges of the lightest neutralino 
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FIGURE 1. (a) Scatter plot of the X\ singlino composition (I./V15I 2 ) versus %\ mass (M^o) for various models. The crosses 

represent the Z'-ino composition (lA^el 2 ) for the UMSSM. The direct constraints are imposed. The upper (lower) UMSSM singlino 
band corresponds approximately to moderate (large) value of s, with the lightest neutralino being MSSM-like for large s. (b) The 
allowed mass range of %9 after applying direct constraints (model parameters, gaugino mass unification My = M\~ 0.5M2, M v ±, 

ATz) and indirect constraints ((g— 2)^, Sl^oh 2 , Mz', domain wall). The M^o bounds in the nMSSM are intended to be illustrative 
and are not necessarily quantitatively precise. 



els where the physical Higgses exist as mixtures of dou- 
blets and singlets [12]. 

We choose a phase convention in which and the 
VEVs are real and positive. K and the gaugino masses 
can in principle be complex, but we restrict our con- 
siderations to real values. We scan M 2 , fl — 50 ~ 500 
GeV with a step size of 1 GeV 4 , s = 50 - 2000 GeV 
with a step size of 5 GeV, tan/3 = 0.5, 1, 1.5, 2, 10, 50, 
and |k| =0.1 ~ 0.75 with a step size of 0.05 (for the 
NMSSM). The gaugino mass unification relation My = 

Mi = |^jM 2 ~ 0.5M 2 is assumed. 

The bounds that we obtain on the the lightest neu- 
tralino mass are 

53 GeV < M y0 < 248 GeV [MSSM] 
OGeV < Mn < 248 GeV [NMSSM] 

Z\ 

OGeV < Mn < 83 GeV [nMSSM, S-model] 

X\ 

OGeV < Mn < 248 GeV [UMSSM] 

Z\ 

These results are only for positive gaugino masses, but 
negative gaugino masses do not change these ranges sig- 



nificantly 5 . The mass bound dependence on tan /3 is quite 
sensitive in some models. For example, at tan/3 w 1, the 
NMSSM and the UMSSM have massless %\ while the 
nMSSM has an upper M o bound [1]; the MSSM vio- 

lates the LEP2 m/, constraint at this tan j3 . The condition 
for the massless (or very light) neutralinos can be eas- 
ily satisfied in extended MSSM models without the fine- 
tuning which is required for the MSSM case [1]. 

Figure 1(a) shows the M o dependence on the singlino 

composition of X\, Wis\ 2 (and also the Z'-ino composi- 
tion |iVi6| 2 for the UMSSM). The dashed lines are the 
MSSM bounds. Singlino dominance is typical in the 
% x of the extended MSSM models. Especially, when 
the Mn is much smaller than the MSSM lower limit 6 
(M o ~ 50 GeV), the singlino is always the dominant 
component. 



5 With Mi = -50 ~ -500 GeV, the %{ mass ranges are M^o = 39 ~ 

254 GeV (MSSM), M % a = ~ 254 GeV (NMSSM), M % a =0.4 ~ 96 

GeV (nMSSM), M y0 = 39 ~ 254 GeV (UMSSM). 
, *■ 
The relaxation of the gaugino mass unification requirement would 

allow the MSSM also to have a very light (bino-dominated) neutralino. 

For the supernova constraint in this case, see Ref. [13]. 



4 We exclude very small Mi or \x. to avoid very light non-singlino 
states. These are also excluded by the chargino mass constraints. 



DISCUSSION 

In the nMSSM (also in the S-model in the decoupling 
limit), because of the small X\ mass, the Z pole is the 
dominant channel [11, 14]. To obtain an acceptable relic 
density with only the Z pole annihilation contribution, 
the lower M^o bound is M v o > 30 GeV which is allowed 

for only small tan/3, while the solution for a 2.4c devi- 
ation from the SM of (g — 2) M favors large tan/3, which 
allows only small M^o. Nonetheless, a common solution 
was found to exist in this model [15]. We refer to Ref. 
[16] for the interesting physics associated with the light 
pseudoscalar Higgs boson including an alternative relic 
density annihilation channel. 

The NMSSM is disfavored by the non-observation 
of cosmological domain walls predicted by the discrete 
symmetries of the model. However, there is an approach 
to interpret the domain wall as the dark energy [17]. 

In the UMSSM, the mass of the new gauge boson Z' 
is also an important constraint. The Higgs singlet VEV 
s should be large enough to satisfy the experimental M z i 
bound, although there are some ways to get around it. 
(See Ref. [1] for the discussion.) 

Approximate lightest neutralino mass ranges in the 
models considered are illustrated in Figure 1(b). The 
dashed regions are disfavored by the indirect constraints 
((g - 2)u, £l r oh 2 , M z i, domain wall). After the M z i lower 
bound is imposed, the UMSSM bound becomes similar 
to the MSSM. In the case of the nMSSM there is a 
tension between the (g — 2) M constraint which favors 
small M o (or large tan/3) and the relic density constraint 

which favors large Mo (or small tan/3). 

CONCLUSION 

Although Supersymmetry at the TeV scale is well- 
motivated, the MSSM is just one of the possible real- 
izations. In fact, the fx -problem suggests that the MSSM 
is incomplete. The solution to the fi -problem suggests 
that an appropriate direction to extend the MSSM is to 
have an extra Higgs singlet whose VEV gives the effec- 
tive /i-term of the EW scale. Extensions of the MSSM 
have extra neutralinos, and the composition of the light- 
est neutralino involves extra components beyond those of 
the MSSM. Because of this, both the mass and couplings 
of the lightest neutralino are modified from the MSSM. 
The lightest neutralino (jj 1 ) is interesting both in particle 
physics (as the LSP) and cosmology (as the CDM), and 
it is therefore important to study and compare properties 
of the Xi in extended MSSM models. 

The properties of the CDM particle, even if it is 
the lightest neutralino, may be quite different from the 
MSSM prediction. For instance, it could be extremely 



light and/or dominated by the singlino, which does not 
directly couple to SM particles except Higgs doublets. 
Similar distinctions of models may occur in the Higgs 
sector. Even if a low energy Supersymmetry is correct, 
its realization may depend on the model. The measure- 
ment of the mass of the lightest neutralino and the de- 
termination of its couplings will be particularly useful in 
testing the MSSM and its extensions at colliders [18]. 
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